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ABSTRACT 



q \ We discuss the role of ambipolar diffusion for simple reconnection in 



a partially ionized gas, following the reconnection geometry of Parker and 
Sweet. When the recombination time is short the mobility and reconnection 
of the magnetic field is substantially enhanced as matter escapes from the 
reconnection region via ambipolar diffusion. Our analysis shows that in the 
t**« ' interstellar medium it is the recombination rate that usually limits the rate of 



reconnection. Consequently, the typical reconnection velocity in interstellar 



medium is ~ (^/T reC omfe) 1 ^ 2 (l + x~ l (3~ l ) : where r\ is the ohmic resistivity, 
x~ x is the ionization fraction, and (3 is the ratio of gas pressure to magnetic 
pressure. We show that heating effects can reduce this speed by increasing 
the recombination time and raising the local ion pressure. In the colder parts 
of the ISM, when temperatures are ~ 10 2 K or less, we obtain a significant 



O 

enhancement over the usual Sweet-Parker rate, but only in dense molecular 

clouds will the reconnection velocity exceed 10~ 3 times the Alfven speed. The 

ratio of the ion orbital radius to the reconnection layer thickness is typically 



a few percent, except in dense molecular clouds where it can approach unity. 
We briefly discuss pro; 
astrophysical plasmas. 



x 

We briefly discuss prospects for obtaining much faster reconnection speeds in 



Subject headings: Magnetic fields; Galaxies: magnetic fields, ISM: molecular 
clouds, magnetic fields 



1. Introduction 

Understanding the mobility and reconnection of magnetic fields in a conducting 
medium is critical to understanding the evolution and origin of large scale magnetic fields in 
astrophysical objects such as stars, galaxies, and accretion disks. Standard dynamo theory 
relies on the turbulent transport of magnetic flux to move the field lines and, implicitly, 
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to change the large scale topology of the magnetic field. However, the substitution of 
turbulent transport for microscopic diffusion is difficult to justify on theoretical grounds. 
Normally we would expect flux freezing to be an excellent approximation to the motion 
of the magnetic field in a highly conducting fluid. Getting field lines to pass through one 
another, or rearrange their connections, must ultimately involve ohmic diffusion. Moreover, 
as the field lines are stretched, they will strengthen and exert forces that will prevent further 
stretching or deformation. Recently several different authors (e.g. Cattaneo and Vainshtein 
1991, Gruzinov and Diamond 1994) have argued that this raises an insuperable obstacle to 
turbulent transport, and standard mean-field dynamo theory, in highly conducting fluids. 
One of us QVishniac 1995|) has proposed that this problem may be solved through the 



formation of intense flux tubes, but this assumes that reconnection is rapid enough to allow 
the formation of such structures in a small number of dynamical time scales and that the 
plasma has a high f3 so that the magnetic field can be distributed intermittently. 



Estimates of reconnection speeds based on a simplified geometry ( jSweet 19581 , Parker 
19571) gi ve 

V rec *(^) 1/2 = v A Re^ 2 (1) 

where L x is the typical structure scale, va is the Alfven speed, r\ is the ohmic diffusion 
constant, and Res = (vaL x /tJ) is the magnetic Reynolds number. In general, we expect 
that va is comparable to the local turbulent velocity, so that this speed will be many orders 
of magnetic slower than typical fluid velocities in astrophysical objects, where Res >> 1. 
However, we note that observations of magnetic fields in the solar chromosphere and corona 
(cf. Dere 1996, Innes et al. 1997) suggest that reconnection can occur at speeds ~ O.Iva- It 
is not clear what conditions are necessary for such rapid reconnection or what mechanism 
is operating that allows it. 

In some ways, these problems are especially severe when we consider the evolution of 
the galactic magnetic field ( [Kulsrud fc Anderson 1992] ). The huge scales involved limit 
the number of dynamical time scales available for generating a large scale field and strong 
theoretical objections have been raised against the possibility of a primordial field strong 
enough to eliminate the necessity for a galactic dynamo. In addition, the magnetic Reynolds 
number in the galactic disk is of order 10 20 , so that the Sweet-Parker reconnection speed is 
negligiblely small. 

On the other hand, the interstellar medium is only partially ionized, and in dense, 
cool clouds the ionization fraction is very small. This raises the possibility that magnetic 
field lines which move with the ion velocity may nevertheless be capable of a high degree 
of mobility relative to the bulk of the gas. Ultimately, the speed of magnetic reconnection 
is limited by the width of the current sheet dividing magnetic fields of sharply different 
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orientation. Although importance of ambipolar diffusion has been long stressed (see 
Mestel 1985), the results obtained so far are contradictory. Recent one dimensional 
numerical work by Brandenburg and Zweibel (1995) and Zweibel and Brandenburg (1997) 
suggests that ambipolar diffusion can give rise to narrow current sheets, whose widths are 
ultimately determined by the size of particle orbits in the plasma. This in turn suggests 
that reconnection in the neutral parts of the interstellar medium may play a critical role 
in the galactic dynamo (cf. Subramanian 1998). At the same time analytical studies of 
reconnection in the presence of ambipolar diffusion by Dorman (1996), also based on a one 
dimensional model, led to a different conclusion. 

These ambiguities motivate our current study. In this paper we reexamine the role 
of ambipolar diffusion in Sweet-Parker reconnection, including the effects of ion pressure 
and finite recombination rates for the ions, and allowing for the transverse loss of neutral 
particles during reconnection. We calculate reconnection velocities and show that for most 
astrophysically important cases reconnection is limited by the recombination rate. We 
find a substantial enhancement of the reconnection rates even in media with low rates of 
recombination. This enhancement is not enough, however, to account for magnetic flux 
tube formation in the ISM (as discussed in Lazarian & Vishniac (1996) and Subramanian 
(1998)). At the same time the rates obtained here can explain topological changes that 
accompany stellar formation, that is, the disconnecting of the magnetic field of a collapsing 
cloud from the interstellar magnetic field. 

In section II of this paper we will briefly review the physical basis for the Sweet-Parker 
rate, and for the enhancement due to ambipolar diffusion. In section III we estimate the 
effects of ion pressure and give a simple mathematical model supporting our claims. In 
section IV we calculate rates of reconnection for different phases of the interstellar medium. 
Finally, in section V we discuss the implications and limitations of our work, and summarize 
our conclusions. 



2. Sweet-Parker Reconnection and Ambipolar Diffusion 

2.1. The Case of the Imperfect, Ionized Gas 

We start by reviewing the physical basis for the Sweet-Parker reconnection rate. This 
is well known material, but including it here will make subsequent arguments concerning 
ambipolar diffusion and recombination clearer. The simple geometry that forms the basis 
of the Sweet-Parker reconnection rate consists of two regions of opposite magnetic polarity 
facing one another. The y axis is perpendicular to the field lines and gas is approaches the 
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midplane with a velocity — u y for y > and u y for y < 0. The magnetic field lines curve 
apart a distance ±L X along the x axis. As the field lines curve apart the gas streams to 
the right and left with a velocity u x . The magnetic field is zero along the x axis, but rises 
to ±B a distance L y on either side. This geometry is illustrated in figure 1, although the 
arrows for gas flow are appropriate for the case of ambipolar diffusion. When the gas is 
completely ionized the flow along the x axis is confined to the layer of width A. 

The evident neglect of the existence of a third dimension is the major weakness of this 
picture. However, the tension of the magnetic field lines, and their consequent resistance to 
bending, makes the role of this neglected dimension quite complicated. Here we will simply 
assume that this geometry provides a useful approach to the problem of reconnection, and 
defer consideration of more complicated geometries to a subsequent paper. It is important 
to keep in mind that in three dimensions the magnetic field along the reconnection surface 
does not actually vanish, but simply reduces to a component in the z direction which is 
common to both magnetized regions. This does not effect our arguments in any way, but 
implies that one should resist the temptation to think of the plasma in the reconnection 
zone as essentially unmagnetized. 

Conservation arguments can be used to estimate u x and constrain the reconnection 
geometry. First, we note that the pressure along a line parallel to the y axis and passing 
through the origin is approximately constant, i.e. 

B 2 

P H m constant. (2) 

8ir 

At the midplane the magnetic pressure vanishes. Following the fluid to the right or left, one 
comes (in a distance ~ L x ) to a region unaffected by the magnetic pressure. This implies 
a pressure excess in the reconnection region of B 2 /8tt. This pressure excess is sometimes 
derived by considering the dissipation of magnetic energy in the reconnection layer, but this 
gives rise to the mistaken notion that an efficiently radiating reconnection layer will not 
have a gas pressure excess. From Bernoulli's theorem 

p 

ul H ~ constant , (3) 

p 

which implies that if u x — at the origin then at x — ±L X 
where va is the Alfven velocity of the magnetized regions. 
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The mass influx into the reconnection region is just 2pL x u y , while the mass efflux is 
2pL x u y = 2pL y u x , so from conservation of mass we have 

Uy = -j^VA . (5) 

We note that the ions need not trace the movement of the magnetic flux in an imperfect 
fluid. The magnetic field current is 



'B 



B-r/VB , (6) 



where the subscript i refers to the ions. Given the geometry used here, the effective velocity 
of the magnetic field can deviate from u y by an amount f]/ L y . To put it another way, the 
magnetic field is no longer flux frozen if L y is sufficiently small and regions of opposite 
polarity can annihilate much faster than matter can be dragged into the reconnection 
region. This acts as a kind of regulator for the size of the reconnection region. If it is too 
large, then reconnection slows to a crawl and whatever external forces are pushing the 
magnetized regions together will continue to do so[|. If it is too small, then reconnection 
runs away and broadens the neutral zone. We can therefore take 77 ~ L y u y and obtain 

(VVa\ 1/2 ,_x 

This is the Sweet-Parker reconnection rate given in equation ([]]). 

One important aspect of this model is the existence of a pressure deficit, of order 
the magnetic pressure, at a distance ±L X along the x axis. In a turbulent medium with 
equipartition between the magnetic and turbulent energies such pressure excesses will come 
and go on a time scale not much longer than L x /va, so that this stationary model can be 
no more than a very approximate guide to the structure of the reconnection layer. However, 
the existence of magnetic tension guarantees that if the magnetic field reaches a locally 
persistent equilibrium, then we can naturally expect such pressure fluctuations to last as 
long as the magnetic structure itself. In other words, the tension in magnetic 'knots' is 
sufficient to guarantee a downstream pressure deficit of order B 2 /8ti. One might suppose 
that the accumulation of ejected ions will erase this pressure deficit, but as reconnection 
proceeds the reconnected magnetic field lines will be pulled away from the reconnection 
region, bearing with them the ejected ions. It is also true, but not obvious, that ejected 
neutrals can usually be ignored in this model, a point we will return to later. 



lr The role of the external forcing is, in fact, more complex. It can squeeze the material from the 
reconnection zone and thus increase the local Alfven velocity. 
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2.2. The Case of the Imperfect, Partially Ionized Gas 

When the fluid is partly neutral the charged and neutral particles will no longer move 
together, and Vj is no longer the same as the bulk velocity of the fluid. We can estimate 
Vj — v n by balancing the pressure exerted on the ions with the collisional drag due to 
collisions with the neutrals. We get 

V5 2 2 V n -Vi V n -Vj 2 

— — + c^pi = pi— = p n — = -c n Vp n , (8) 

where ij >n and t ni are the collision times for ions impacting on neutrals and vice versa, 
respectively, and Cj and c n are the ion and neutral sound speeds. If we neglect the ion 
pressure gradient (more on this later) we can write the neutral particle velocity as 

V£? 2 

V n = Vi + t n>i , (9) 

8np n 

or 

VB 2 , N 

V* = V - -T, tn,i , (10) 

8np 
where p and v are the total density and bulk velocity of the gas. This implicitly assumes 
that the neutrals are moving relatively slowly, and dominate the fluid mix, so that the 
neutral pressure gradient can balance the magnetic pressure when the ion pressure gradient 
fails to do so. We can use this expression in equation ([]) to write the magnetic field current 
as 

J B = v(g)B-^ + ^yVB , (11) 

where va is the Alfven velocity relative to the total gas density. From here on we will 
define r\ ara u = v 2 A t n ^. Typically va ~ 10 5 cm/sec. Since t n> i ~ 4.8 x 10 8 /rii seconds at low 
temperatures ( Palgarno 19581 ), an d T] is of order 10 9 cm 2 /sec, r] ambi will be many orders of 
magnitude larger than rj in cool, low density regions in the interstellar medium. 

Equation (|TT|) seems to imply that we can use r\ amhi 3> r\ in place of the usual 
expression for Ohmic diffusion. The neutral gas will stream outward in the x direction with 
a velocity ~ va, driven by the pressure gradient created by collisions with the in-flowing 
ions. In reality this derivation includes assumptions which impose severe constraints on 
the substitution of r/ am bj for rj. The most important is that we have assumed that the ion 
pressure gradient is negligible compared to the magnetic pressure gradient. However, if the 
ions move inward at a steadily increasing speed towards the reconnection surface then for 
y <^ L y the only escape route for the accumulating particles is to join the opposing flow of 
neutrals. This can only be true if the recombination time is short, so that ions and neutrals 
change identity easily. This expression is commonly used in calculating the ambipolar 
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diffusion time for molecular clouds, in which case one only requires the recombination time 
to be no greater than the time it takes the magnetic field to drift out of a cloud. 

Here we are concerned with the formation of a reconnection layer, in which the relevant 
time scales are much shorter. We begin by assuming rapid recombination. In that case 
we can see that the width of the reconnection region, L y is set by the matter diffusion, as 
above, by equating the drift velocity to u y . This implies 

"•~(-etj = k~tt) • (12) 

when rjambi > 77. Later on we will refer to this expression for u y as V max to emphasize that 
this is the maximal velocity of reconnection that is obtainable through ambipolar diffusion. 
By the same reasoning we have 

1 /9 

Ly ~ f ^ mbiL A . (13) 

In obtaining expressions QT2D and flT5| ) we assumed, first, that the outflow from 
the region L y happens with a velocity va and, second, that on leaving the reconnection 
layer neutrals diffuse slowly out while moving together with ions over the distance 
comparable with L x . The former assumption depends on the pressure differential around 
the reconnection region, which is assumed as part of our basic geometry, and the dynamics 
of the outflow, which we check below by studying the structure of the ambipolar diffusion 
layer. The latter assumption constrains the geometries of the reconnection layers considered. 

It's important to note that the actual process of reconnection takes place only in 
a narrow zone where rj > r] am bi- However, unlike the usual Sweet-Parker result, the 
width of this inner zone does not determine the reconnection speed. Instead its width 
automatically adjusts itself to match the reconnection speed determined by ambipolar 
diffusion. Ambipolar diffusion removes matter from the reconnection zone and enhances 
the reconnection velocity. 

The existence of this broad outflow raises the possibility of a problem in our model. 
Since the gas outflow is not confined to the layer where the magnetic field is actually 
recombining, we cannot assume that the ejected gas will be removed from the downstream 
flow by the relaxation of the reconnected field lines. Instead, the neutral component of the 
gas has to diffuse into the wedge of reconnected field lines by overcoming the ion-neutral 
drag. Since the downstream pressure deficit will help push the neutrals in this direction, 
this diffusion will occur at a rate ~ f] a mbi/L y . In the mean time the transverse motion of the 
neutrals will be approximately along the field lines, since the ions and neutrals are strongly 



coupled through collisions and move together at a speed ~ va- Consequently, the gas will 
move a distance VALy/rj ambi ~ L x while diffusing onto the reconnected field lines. It is a 
basic condition of this model that the reconnection region is embedded in some larger flow, 
allowing room for the relaxation of field lines and the escape of ejected ions. Since this 
distance is no larger than the size of the reconnection region it is a short enough that we 
can assume that the accumulation of ejected gas does not pose a problem for the continued 
flow of gas from the reconnection region. 

We can estimate the width of the resistive zone by considering the motion of the field 
lines for y < L y . Within this zone the magnetic field lines speed up as the pressure gradient 
steepens. Since the bulk velocity remains comparable to V rec , while the magnetic field 
strength plummets, the magnetic flux is carried by the second term on the LHS of equation 
(|TT|). Consequently, we have 

for y <C L y and where B^ is the magnetic field strength far outside the reconnection layer. 
This has the solution 

B = B 00 (^^y + c) 13 , (15) 

where i] am ^ i00 is the ambipolar diffusion coefficient far from the reconnection zone and C is 
a constant equal to B(y = + ) 3 /i?^ o . (The magnetic field strength near y = 0, but outside 
the layer where reconnection actually occurs, is defined as B(0 + ).) 

The ion velocity, which is also the inward speed of the magnetic field lines, is given by 

/ 3y rec X" 1 / 3 

Vi = Vambidy In B = V rec r -^-y + C\ . (16) 

\']ambi,oo / 

If C = then this expression diverges near y = 0. One resolution of this problem is to note 
that Vi is limited by the local value of va/x 1 ^ 2 , where x is the ionization fraction of the 
gas. As Vi approaches this value the assumption that the acceleration term is negligible is 
violated and we have Vi — > va/x 1 ■ This implies a limit on C, which can be used as an 
estimate for B(y = + ). Using equations ( |T5| ) and ( JIB]) we find 

B(0 + ) ~ Boo ( ^^] =B 00 ( XVA T ootn > t ) 1/ ' . (17) 

This expression will be relevant provided that rj remains less than 7] am bi as the magnetic 
field approaches this asymptotic limit. More precisely, this is relevant when 

/ \ 1/2 

V < Vambi,oo — T • (18) 

\ J^x v A,oo / 
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In this case the asymptotic value of Vi is 

^(0 + )=^ iOO (^^) 1/4 , (19) 

and the width of the resistive reconnection region is 

A«-^ T «fMy /2 x i/4f_?_y / v_5_y /4 . (20) 

Vi(0+) \VA,ooJ \Vambi,ooJ \L x V A ,ooJ 

In either case, we expect that B(0+) <C B^ and consequently that magnetic pressure 
scales as y 2//3 throughout the outer layers of the reconnection region. This implies that the 
overpressure in the layer L y scales as -8^,(1 — (y/L y ) 2 ^ 3 ) and that the outflow velocity will 
be of the order of va- This justifies a posteriori our assumption of rapid outflow in this case. 



When the condition expressed in equation ( |T8| ) is not satisfied V{ stays below the ion 
Alfven speed at all times. In this case we can take C = and determine the width of the 
resistive reconnection region from the condition that 



^«fo(l/ = ±A)| . (21) 



From equation (jI5[) this implies that 



A^(^L] 1/2 ^L- . (22) 

\vA,oo J "Iambi, oo 



We see from equations ( BO ) and ( B2] ) that the width of the resistive reconnection region 
is much smaller in this case than when ambipolar diffusion is negligible. This is expected, 
since a narrow reconnection region is necessary for a fast reconnection speed, but it also 
raises the question of local heating. The bulk of the magnetic energy is dissipated through 
expansion as the ions speed up, there is an irreducible minimum which is annihilated inside 
the reconnection region proper. The volume heating rate will be 



£^B(y = + r^L ( (23 ) 



but since A = r)/vi this is 



g&B(y= + ) 2 ^ = B 2 J^ = B 2 J^ Vambl ' QO . (24) 

T) 1] L x 1] 

We conclude that the local heating rate is significantly enhanced in the center of 
the reconnection region, in the resistive layer, and that this may introduce difficulties for 
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a model which relies on an abundance of neutral particles to carry mass away from the 
reconnection region. We will not attempt to calculate the consequences of this heating 
here. It is probably not the most severe constraint on ambipolar diffusion in the interstellar 
medium. Instead, we suggest that the long recombination time for ions, coupled to the 
effects of ion pressure, poses a much greater problem. We examine this point in the next 
section of this paper. 



3. The Role of Ion Pressure 

The critical assumption in the preceding section is that it is reasonable to assume 
that the ions and neutral particles maintain their usual ratio within the reconnection layer. 
We can see why this important by restricting our attention to the resistive layer. The net 
inward flux of ions will be > V rec n itOQ . This constitutes a lower limit, since within the whole 
reconnection region of width L y the magnetic field lines speed up, so that their density 
drops, and the same effect will depress rii below its equilibrium value. A uniform neutral 
distribution will therefore add incoming ions to this flow. 

The rate at which ions leave the resistive region through recombination is 

^ (0) 7n\ =n ^ \,(n W 77^> n ^Vrec ■ (25) 

'recombyyj vi\y+)'recomb\y) 

This assumes that the loss of ions through their expulsion in the x direction is 
negligible, but if that process dominates then we will recover the usual Sweet-Parker result. 
When the recombination rate is small, ions will accumulate in the resistive layer creating 
a strong ion pressure gradient, and the reconnection process will be limited by the rate at 
which ions can recombine with electrons within the resistive layer. 

We can understand the role of the ion pressure by equating the ion density flux with 
the recombination losses in the resistive layer. We have 

Vi(0 + )Pi(0 + ) = -^£r ■ (26) 



<recomb\y ) 



This implies that 



V2 /rrm at/ B 2 \ 1/2 



M0+) = (-mj\ (m* + ^^^ i , 2 7) 



\Trecomb(0)J \ T(0) ^(0+) 87Tp i)OO C?(0) , 



where we have invoked pressure balance for the ion density in the resistive layer, and 
assumed that B/pi is a constant outside of the resistive layer. 
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The recombination rate is an e} where n e ~ rij and a is the recombination coefficient 
(see Spitzer 1978, p. 107), which scales approximately as 2 x 10~ 11 </>(T)/T 1 / 2 cm 3 s _1 for a 
hydrogen plasma, where is a slowly varying function of temperature. In an ionized plasma 
the resistivity is ~ 2.65 x 10 12 T~ 3//2 , and this expression should remain valid as long as the 
scattering length for electrons is determined by inelastic collisions with ions. In a largely 
neutral gas the exact dependence on density and temperature can be fairly complicated. 
However, for our purposes a crude estimate is all that is required and in any case the 
density of ions in the resistive layer is greatly enhanced in order to balance the magnetic 
pressure in the surrounding plasma. Using these results we can rewrite equation ( p7]) as 

/ 2,(0+) \ " 2 (no+n 3 ' 2 1 v„ bi \ 

A + '~'W mk (0 + )j \T(0) ) {' + t.,(0 + ) 8^,^(0+) ) ' [Z01 

Equation (^) is only an upper limit on V rec , which in principle could be much lower if 
there exists an outer layer dominated by ion-neutral drag. However, we can show that this 
is unlikely to be the case. If such a layer exists then we can see from equation ( Jl6| ) that the 
ion velocity in this layer becomes 

\Vambi,oc \Vi{0 + )J J 

This in turn implies an outer layer width of 

T r lambi,ooV rec for\\ 



Vi(0- 



and a reconnection velocity of 



V rec =(^^\rM^ . ,31, 



J x 
\ vA'lambi,c 



This could represent an enormous reduction of V rec from its upper limit of t>j(0 + ), but this 
outer layer will not form unless the cumulative drag from the neutrals in this layer is large 
enough to lead to a significant reduction in the local magnetic field strength. Otherwise the 
outer layer will collapse and fj(0 + ) — > V re 



rec- 



At large distances the ions and neutrals move together with a velocity V rec towards 
the reconnection zone. The accumulation of neutrals will lead them to decelerate and 
accumulate in a layer of width L n . If we assume that this layer is much thicker than the 
actual zone of reconnection, then the relative ion-neutral velocity in this layer will be ~ V r 
and the net drag force will be 

PnVrec—^ ~ 8P n , (32) 



rec 
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where we have assumed that the drag force is balanced by the excess neutral pressure. The 
neutrals will be expelled laterally at a speed v x ^ n given by 

V x,n ~ -— ■ (33) 

Pn 

Using Eqs. (0), (0) and conservation of mass we conclude that 

/ V rec lj x \ 



(34) 
and 



m,i 



Ln-M^^y 73 - (35) 



J x 



The condition that the total drag within this layer has a negligible effect on the reconnection 
layer is SP n /p n <C v\ or t> xn <C va (which is also the condition that the neutral ejection 
velocity be much less than va)- This is equivalent to 

\r2 „ Vambi v A /„ R \ 

Kec < — ; • (36) 

i-'x 

In other words, as long as the reconnection speed is less than we would obtain from the 
naive substitution of r\ amhi for r\ in the usual Sweet-Parker formula, we can ignore the 
ion-neutral drag outside the reconnection layer and set Vi(Q + ) equal to V rec in equation (p 



A significant complication is that the momentum of the ejected neutrals may diffuse 
over a distance greater than the value of L n given in equation (|35l), thereby increasing L n 
and possibly the role of the neutral drag in the reconnection layer. This effect will become 
important when 



- (37) 



L n L x 

where v n ~ c 2 s t ntn is the neutral gas viscosity. In this limit we need to replace equation ( p3|) 
with 

"mi Vn /oo\ 

— T PnJjV X:n . (38) 

L x J^n 

If we combine this result with equation (R2f) and the condition that V rec L x ~ v x n L n we get 



L n ~(cXn*n,i) 1/4 4 /2 - (39) 



The condition that the ion-neutral drag does not dissipate a significant fraction of the 
magnetic energy outside the reconnection layer is v\ 3> V rec L n /t nj i, which can be rewritten, 
with the aid of equation (|3UJ) as 

/r2„2 f \V4 
V\ » Vrec \^f^\ , (40) 
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or, using the definitions of V max and r] ambi , 

Kec«K laa; (^) 1/2 f^y 4 . (41) 

For largely neutral gases the ratio of t ni to t n ^ n will be of order x -1 , and va is usually of 
order c s . Consequently, under typical conditions in the interstellar medium equation ([41]) 
will be satisfied by a comfortable margin (cf. Table 1). 

Once again we need to consider whether or not the broad, slow outflow of neutrals from 

the reconnection region will result in an accumulation of gas outside the reconnection region, 

along magnetic field lines that have not yet undergone reconnection. The appropriate 

diffusion coefficient is, as before, v A t n ,i- The diffusion scale is L n , the thickness of the 

stagnation layer in the neutral flow. We see from equations (|i~2"D, (|34"D, and (^) that 

this implies that the transverse distance covered by the neutrals while diffusing into the 

reconnected wedge will be 

L 2 fV\ 4/3 

W^U^) • (42) 

Since this is always less than L x we conclude that the ejected neutrals can be ignored. 

This case is illustrated schematically in figure 1, with the local velocities of the neutrals 
and ions indicated by thin and thick arrows. The material in the reconnection layer is 
actually moving in the x direction with a speed ~ va, but since most of the gas leaves the 
reconnection layer in the y direction as neutral particles we have ignored the motion of the 
ions in that layer. 

Our only remaining worry is that the transverse expulsion of neutrals may serve to 
remove large numbers of ions from the reconnection zone. It is certainly reasonable to 
assume a tight coupling between the ion and neutral transverse speeds, since the transverse 
shear rate, v x>n /L x is much smaller than the ion coupling rate, (p n /Pi)tn}- However, for 
typical ISM parameters we also have T~^ comh > va/L x > v x>n /L x , so the ion fraction in the 
gas is maintained even while individual ions are expelled. We also note that even if we 
ignored recombination the equations of continuity for ions and neutrals combine to give 

Vid y In pi - v n d y \np n + V(^ - v n ) = 0. (43) 

Within the outer layer of thickness L n the neutral particles come to a halt and V{ remains 
close to V rec . Consequently the density scale height for the ions is not less than L n and 
some large fraction of the ions will reach the resistive layer, despite transverse losses. 

We conclude that 

K _(^n^n 1+ ii , (44) 
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where (3 is the ratio of the gas pressure to the magnetic pressure. The width of the 
reconnection layer is 

A = ^ = (r /oo T recomb , 00 )Wl + ^) . (45) 

* rec \ PX J 

The principle source of ion heating in the reconnection layer is ohmic heating due to 
the dissipation of the magnetic field. Cooling can take place either through radiative losses 
or through the transfer of energy to the neutral gas. The neutrals will, in turn, shed their 
excess thermal energy through radiative losses and/or by escaping from the reconnection 
layer and diluting the extra heat over a large volume. Since neutral mean free path is 
typically much larger than A, and since c s ^> A/t n ^, the ion cooling rate will be controlled 
by the ion-neutral collision rate. In other words, 

where we have assumed that T(0) 3> Too, since otherwise the correction is of no interest, and 
used Boo since we have already shown that B^ ~ -8(0+) for reconnection in the interstellar 
medium. We have also assumed that the electrons and ion share the same temperature 
in the reconnection layer. Given that the charged particle pressure in the ionization zone 
balances the external magnetic pressure, equation (fill) implies 



V 2 3 

rec ^ 



77(0) t it n 

We can rewrite this with the aid of equation (03) as 



(47) 



/ 3r recomb]00 \ 0(Tqq) / T(0) \ 
W ^ t hn ) 0(T(O)) [t^ ) ■ [ j 

This implies that when reconnection layer heating is important we should substitute the 

expression 

Crn \ 3/2 / 1 \ -6/5 

^y) « 31 ^1 + ^J x'^(T(0)r 3/5 T°J, (49) 

into equation (H3). We see at once that T(0) > T^, so that this correction is appropriate, 




only in cold regions with x ^C 10 3 (assuming (3 is of order unity). As we reach this limit 
we go from V rec oc x -05 to V rec oc x 01 . 

It is important to remember that for strictly one dimensional formulations of this 
problem there is no natural choice for L n and the actual reconnection speed will remain 
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sensitive to the computational box size and/or the boundary conditions. Although L x does 
not appear in equation ( f44|) it is still present implicitly as part of the constraints that make 
this a well-posed problem. Purely one dimensional calculations will not capture all the 
relevant physics, and may result in a wide variety of estimates for V rec and A. 



Equation fl44|) constitutes the main result of this paper, and may be regarded as a 
generalization of the Sweet-Parker reconnection rate to an ambipolar medium with a long 
recombination time. The major effect we have neglected here is additional ionization within 
the reconnection layer caused by the reconnection process. We have included the effect of 
local heating on the resistivity, gas density, and recombination rate. We note that even 
taking T(0) = T^ we can show that reconnection in this model is rather slow. 



We should note that equation (|j) is correct only if r recom b is much less than the 
shearing time L x /va- In the opposite limit, ions will escape from the sides of reconnection 
zone rather than through recombination. In this case we are almost back in the regime 
described in section 2.1. However, there are two important physical effects which have to 
be considered. First, we have already seen that the partial decoupling of ions and neutrals 
implies a large concentration of ions in the reconnection layer itself. Consequently the rate 
of ion ejection is enhanced by the ratio rii(0)/ni }OO . The ion conservation equation becomes 

Av e j ect rii(0) = L x V rec n iy00 . (50) 

Second, the neutrals will diffuse out the reconnection layer, spreading the transverse 
momentum and creating an additional drag on the ions. Since these effects work in opposite 
directions, it's not immediately obvious whether the presence of neutrals increases or 
decreases the reconnection speed when T recom b is very large. 

Balancing the ion pressure gradient in the x direction with the neutral drag, we get 

-r- « Pi(0)-— ■ (51) 

The transfer of momentum to the neutrals is then balanced by the viscous drag on the 
neutrals, so that 

Pi(0)— - — A « pnVejecty-, (52) 

where we have assumed that Av x , the difference between the ion and neutral particle 
velocities in the x direction, is small compared to v e j ect . The width of the spread of x 
momentum is given by equation (|37| ) so we conclude that 

v 2 A A n (v 3 eject L m u n ) 1/2 . (53) 
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Remembering that A ps rj(0)/V rec we can combine equations ( (50]) and (|53|) to obtain 



^(M^fV^y'v^r ( 54) 

\ L x ) \ v n ) \ ra ij00 / 
Consequently, we can estimate the corrected Sweet-Parker reconnection velocity as 

1/2 /„ \l/8 / T \ 21 / 16 / 1 \ 3 / 8 



Vn 



mm m ^ ■ 



In spite of the fact that typically z/ n ^> r/oo, this estimate will usually be larger than the 
standard SP reconnection speed for largely neutral, magnetized gas. 

Finally, we note that our estimates are based on the assumption that the reconnection 
layer is in a steady state, meaning that the local pressure excess is comparable to the 
magnetic field pressure and that the accumulation of ions in the reconnection layer has 
reached a steady state. If the magnetic field lines are being compressed by turbulent forces 
then the former condition will be reached in about one eddy turnover time, or, assuming 
Va ~ Vturb, about L x /va- Otherwise we should replace the eddy turn over time with the 
large scale time dynamical time scale. At earlier times the confining pressure will be less 
than our estimates, and the reconnection rate slower. On the other hand, filling up the 
reconnection layer with ions will take a time ~ (A/V rec )(ni(0)/ni tOO ) and at earlier times 
the reconnection layer will be compressed to a smaller width, with a consequent increase in 



the reconnection rate. We can rewrite this time scale using equation (44). We obtain 
A ru(0) 7/(0)^(0) 



V rec Tti ! oo ' rec "*»oo 




Trecomb,oc -ffT~ U + ~ « " ( 56 ) 



This will usually be somewhat less than an eddy turnover time in the interstellar medium. 
We conclude that the early phases of reconnection are apt to be characterized by weak 
compression and slow reconnection, and in any case that the reconnection of magnetic flux 
is unlikely to proceed much more rapidly in the early phases of reconnection. 



4. Phases of the ISM and reconnection 

Does ambipolar diffusion actually lead to a dramatic enhancement of reconnection 
rates in the ISM? In Table 1 we give our results for idealized partially ionized phases of 
the ISM. We include, for reference, V max , the reconnection velocity when r recorn b — > 0, 
and Vsp, the Sweet-Parker reconnection speed without corrections for natural drag or ion 
density enhancement. More realistically, we also show V rec , in which the reconnection speed 
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is limited by the recombination rate, and Vcsp, the corrected Sweet-Parker reconnection 
speed. For molecular hydrogen clouds the reconnection layer heating is important in V rec , 
and has been included using equation (^Uj). The reconnection velocities Vsp and Vcsp will 
depend on the assumed values of L x . Here we take L x = 1 pc so that both these speeds 
scale as (1 pc/L x ) 1//2 . We also make use of the evidence that for all the phases shown in 
Table 1 the Alfven velocity Va is of the order of a couple of kilometers per second. Here we 
assume that it is 10 5 cm s" 1 . Finally, each ISM phase we give nf, the mass fraction of that 
phase. A discussion of the values we use here is contained in Draine and Lazarian (1997). 
We note that an additional 22% is contained in ionized gas. 

There are several important points to be noted in connection with Table 1. First for 
all the phases presented in Table 1 it is the "recombination limited" reconnection rate (see 
Eq (PD) that determines the reconnection speed on scales of ~ 1 pc. Apart from diffuse 
warm HI, where reconnection is slow anyhow, the corrected Sweet-Parker reconnection 
speed given by equation fl55|) is at least one order of magnitude slower than V rec . This 
situation changes only at much smaller scales where outflows from reconnection regions can 
carry out most of the ions faster than recombination can eliminate them. 

Second, in all phases of the ISM, V max is clearly a gross exaggeration of the degree to 
which ambipolar diffusion can enhance reconnection. Nevertheless, V rec is very much larger 
than Vcsp f° r a H the cold phases of the ISM. The latter is, in turn, considerably larger 
than Vsp- In both cases, this increase is due to the compression of ions in the reconnection 
layer. The conclusion that V rec > Vcsp is due to the double contribution of the ion density 
enhancement in the former case, that is, not only does it contribute directly by enhancing 
the efficiency of ion loss relative to the incoming flux of ions, but it also contributes by 
raising the recombination rate. This explains why the warm neutral gas shows only a small 
increase in its reconnection speed. It also suggests that weak magnetic fields will have much 
slower reconnection speeds. 

Third, V rec is always much smaller than va- The only case in which it exceeds 10~ 3 x va 
is that of dense molecular gas. Moreover, this phase has a modest mass fraction of about 
3%. While this result implies that reconnection should be efficient enough to remove 
magnetic flux from star forming regions, and disconnect the internal magnetic field from 
its environment, it does not explain how reconnection could be fast enough to sustain the 
galactic dynamo. 

Finally, these calculations are based on the assumption that the plasma can be treated 
as a fluid. Given the low collision rates in the ISM, this assumption requires, at a minimum, 
that all scales of interest be much greater than the ion Larmor radius. In particular, we 
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require that A 3> Rlaou or 

V rec V th,ion?TlionC I * \ / -, , J- \ *th,ion CX tr*-7\ 

VW eB \riooTrecomb,oo ) \ PX ) V A U^ ion 

where w p ^ on is the ion plasma frequency. If we evaluate this expression for the idealized ISM 
phases shown in Table 1 we find that it rises from a low of about 3 x 10~ 2 for diffuse cold 
HI to about 0.13 for diffuse molecular gas. There is a weak dependence on reconnection 
layer heating which increases the ratio for molecular gas. As before, we have calculated the 
local heating using equation P5|). For dense molecular gas the ion Larmor radius is roughly 
the same as the reconnection layer thickness. This suggests that it would be appropriate 
to reconsider this case without the use of the fluid approximation. Such a calculation is 
beyond the scope of this paper. We merely note that our results for this case should be 
treated with caution. However, it seems unlikely that the reconnection rates will decrease 
significantly when we allow for the finite Larmor radius of the ions. If ?%j on ^> A then the 
ions do not feel the magnetic fields in the reconnection layer. Nevertheless, they will be 
confined to this layer by the electric field of the electrons, which are in turn confined by 
the magnetic pressure exerted on the current-carrying component of the plasma. It follows 
that the electrons and ions will be moving together as in the hydrodynamic approximation. 
Possible deviations from this approximation in the form of instabilities, if they take place, 
will only enhance magnetic diffusivity and reconnection. 



5. Discussion and Conclusions 

In this paper we have examined the role of neutral particles on the structure of 
reconnection regions in the ISM. In so doing, we have assumed that the basic geometry 
posited by Sweet and Parker accurately describes reconnection in three dimensions. We 
find that ambipolar diffusion can give rise to much faster reconnection speeds than those 
obtained by the naive use of the usual Sweet-Parker formula, but these speeds are still very 
small compared to the local Alfven speed. In addition, most of the increase is due to the 
compression of ions in the reconnection layer, and would not occur in a fluid where the ion 
and magnetic pressures were comparable. 

How can we interpret these results? Obviously, there are physical effects which we 
have ignored here, but which might affect the geometry or the physics of reconnection. 
Assuming our results are a realistic description of reconnection in the ISM we can examine 
their implications for the structure of the magnetic field in the ISM. We will address the 
latter issue first, and then briefly take up the prospect of obtaining faster reconnection 
speeds under realistic circumstances. 
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What have we learned about flux tube formation? In a previous paper (Vishniac 1995) 
one of us showed that flux tubes should be a natural feature of turbulent high beta plasmas. 
This is consistent with observations of the solar photosphere, but not otherwise directly 
testable. In parts of the interstellar medium the effective plasma beta can be large, if one 
accounts not only for thermal plasma pressure, but for the pressure from compressional 
shocks and cosmic rays (Lazarian & Vishniac 1996). Therefore, it seems possible that flux 
tubes could form in the interstellar medium. The basic mechanism involves the process of 
turbulent pumping, in which flux tubes are stretched by the surrounding turbulence, and 
then twisted to produce close loops of flux which self-destruct and release their entrained 
ions into the surrounding medium. This assumes that reconnection is efficient. 

However, our results indicate that ambipolar diffusion acts against the flux tube 
formation process. In addition to the small reconnection speeds we have obtained, 
ambipolar diffusion makes the tubes leaky since it allows neutrals to cross magnetic field 
lines, diffuse into the flux tubes, and subsequently fill them up with newly produced ions. 
When recombination times are short we can show that flux tubes can still be formed. In 
this case ambipolar diffusivity works in exactly the same way that Ohmic resistivity works 
in a completely ionized plasma. The trade-off between more rapid ambipolar reconnection 
and ambipolar leakage results in a rather weak dependence of turbulent pumping efficiency 
on the ambipolar diffusivity. The situation changes dramatically if the recombination times 
are large, as they are in the ISM. In this case ambipolar reconnection alone cannot provide 
the necessary rates of turbulent pumping. 

We conclude that flux tubes will not form in the ISM unless there are ways to obtain 
dramatically higher reconnection speeds. The only possible exceptions are the dense 
molecular parts of the ISM. On the other hand, we note that the best evidence for the 
existence of flux tubes comes from the solar photosphere, where the high collision rates 
imply that ambipolar diffusion is only slightly more effective than ohmic diffusion. In 
addition, the best evidence for high reconnection speeds comes from parts of the solar 
chromosphere and corona where the neutral fraction is negligible. Evidently there must be 
circumstances where reconnection is much faster than our estimates and these may allow 
the formation of flux tubes in the ISM. 

What do our results say about the galactic dynamo? Our estimates for Vcsp depend 
on L x , and reconnection might occur over small distances at speeds higher than those we 
have given here. On the other hand, we have already assumed L x is of order a parsec, 
and the large scale eddies in the ISM may be as much as two orders of magnitude larger. 
Moreover, it is precisely these large scale motions which are usually thought to drive the 
galactic dynamo. Finally, we note that V rec is insensitive to L x . While it is difficult to 
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be precise about the reconnection speeds necessary to allow turbulent diffusion and the 
existence of a large scale dynamo, the usual assumption is that these speeds must be in 
the range 0.01 to 0.1 times va (see Ruzmaikin, Sokoloff & Shukurov 1988) or the magnetic 
field will become tangled on the smallest scales and show little evidence for large scale 
organization. Our largest estimates for va falls just outside this range and applies to only a 
few percent of the interstellar gas. Once again our results are disappointing in this regard 
and suggest the necessity for some way to obtain much larger reconnection speeds. 

We could consider circumstances where the gas density is much larger, for example, 
a cool stellar atmosphere. The gas in the solar photosphere is largely neutral and dense 
enough that recombination will be very fast. On the other hand, under these conditions 
ti tn is also very much shorter, and ambipolar diffusion is not very much faster than ohmic 
diffusion. For example, if we take n n ~ 10 10 , and B ~ 1 gauss then i] ara u ~ 10 11 , about 10 2 
larger than the ohmic diffusivity in a cool stellar atmosphere. This suggests a reconnection 
speed an order of magnitude larger than the Sweet-Parker estimate, but still quite small 
compared to turbulent diffusivity in a giant stellar atmosphere. 

Are there additional physical effects that might affect our results? We start by 
considering cosmic rays. Cosmic rays are an essential component of the interstellar medium 
and their pressure is roughly of the same order as the pressure of the magnetic field in the 
ISM, suggesting that they may play an important role in all plasma processes in the ISM, 
including reconnection. Our discussion here will follow Longair (1994). 

Cosmic rays diffuse through magnetic field and the magnetic fluctuations on the scale 
of the Larmor radius are the most effective source of scattering. Therefore for our order of 
magnitude estimates we may consider that the time for scattering of cosmic rays is of the 
order t ~ N/ujl, where ul is the Larmor frequency and iV can be defined as the number 
of orbits required to substantially alter the initial motion of the charged particle. In other 
words, if in the course of an individual scattering event the particle deviates over an angle 
$ ~ B rms /B reg , where B rms and B reg are the rms and regular components of the magnetic 
field at the place of scattering, then N satisfies the relation: A rl / 2 $ ~ 1, where the exponent 
1/2 appears above due to the random walk nature of the process. As a result the diffusion 
coefficient for a charged particle with elementary charge e and energy E can be defined as: 

r 2 r cE , , 

t eB 

In the zone of reconnection cosmic rays are streaming along the magnetic field lines and 
perturb magnetic field lines through nonlinear interactions. In these circumstances $ does 
not differ much from unity. As a result values of i] cr can be as large as 10 20 cm 2 s _1 for the 
Mev charged particles. If we assume that the background gas density is negligible, as it 
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may be in the galactic halo, then the Alfven speed approaches that of light, and we obtain 
high reconnection rates ( [Parker 1992| ). 

In the galactic disk, where the density of gas is not negligible, the situation is quite 
different. In general the Larmor radius for cosmic ray particles will be much larger than 
the reconnection layer width. Consequently, cosmic rays cross the reconnection zone in 
a single orbit. Moreover, neutral particles damp the Alfvenic waves, thereby reducing 
the coupling between cosmic ray streaming and perturbations to the local magnetic field. 
Finally, the reduced magnetic field in the reconnection layer decreases the coupling between 
the plasma and the cosmic rays. In the end, the only obvious role for the cosmic rays is 
that they increase the effective pressure outside the reconnection zone, and consequently 
help compress the gas in the reconnection layer. This changes the relation between the 
mass influx and mass efflux so that: 

P\L x u x = p2L y Uy , (59) 

where 

(60) 



P2 _ Pgas 



Pi go-s T" ±cr 

Above Pgas and P cr are the gaseous and cosmic ray pressure respectively. This increases the 
reconnection speed by the square root of the density ratio. Since the cosmic ray pressure 
is roughly the same order as the gas pressure in the interstellar medium, this implies that 
cosmic rays change the reconnection rate just by a factor of order unity. 

How robust is our assumption concerning the basic reconnection geometry? It is well 
known that the kind of narrow current sheet we are positing as the zone of reconnection is 
dynamically unstable to tearing modes flFurth, Killeen, fc Roscnbluth 1963 ). This instability 



persists for reconnection in largely neutral environments, as noted by Brandenburg and 
Zweibel (1995). Turbulent mixing of the current sheets will promote reconnection by 
enhancing the transport of fresh magnetic flux into the reconnection zone and by mixing the 
accumulating ions outward. In fact, at least one author QStrauss 1988|) has proposed that 



this will lead to reconnection speeds of order (5B/B)va and shown that 5B/B will be be 
of order unity in saturated turbulence. Although some enhancement of reconnection speeds 
seems inevitable, there are grounds for doubting that this will dramatically change our 
results. We will defer a detailed discussion until our next paper, but preliminary estimates 
indicate that allowing for the effects of three dimensional nonlinear tearing modes in an 
ionized medium will increase the reconnection speed over the Sweet-Parker estimate by only 
a factor of the magnetic Reynolds number to the one sixth power. 

Finally, and on a more hopeful note, the existence of a third dimension may allow 
the interpenetration of magnetic field lines for favorable geometries, i.e. an interchange 
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instability, which will facilitate reconnection by dramatically increasing the surface area 
of the reconnection layer, (cf. Uchida & Sakurai 1977). We will explore these ideas in a 
subsequent paper. 

We are grateful to E. Zweibel, B. Draine and R. Kulsrud for a series of helpful 
discussions. This work was supported in part by NASA grants NAG5-2858 (AL), 
NAG5-2773 (ETV), and NSF grant AST-9318185 (ETV). ETV is grateful for the hospitality 
of MIT and the CfA during the completion of this work. 
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Table 1. Reconnection speeds for idealized phases of the ISM 

ISM phase r recomh V max V rec V S p V C sp 

diffuse cold HI 4.0 x 10 12 2.3 x 10 3 7.8 9.3 x 10~ 3 1.4 x 10" 1 

n = 30 cm" 3 , x = 10 -3 
T = 100 K, nf « 0.15 



diffuse warm HI 


3.6 x 10 13 


1.9 x 10 3 


4.1 x 10~ 4 


4.3 x 10~ 4 


3.0 x 10~ 4 


n = 0.3 cm -3 , x = 0.15 












T = 6000 K, nf « 0.4 












diffuse H 2 


1.6 x 10 13 


5.6 x 10 3 


(38) 


1.4 x 10~ 2 


0.56 


n = 50 cm -3 , x = 10 -4 












T = 60 K, nf w 0.2 












dense H 2 


4.5 x 10 12 


3.9 x 10 3 


(220) 


3.1 x 10^ 2 


4.0 


n = 10 3 cm -3 , x = 10 -5 












T = 20 K, n/ w 0.03 















Note. - - All quantities are given in cgs units. r recom b is the recombination time outside the 
reconnection layer; Vmax follows from infinitely fast recombination (Eq. (|12|)); V rec allows for 
a realistic recombination rate (Eq. (0)); Vsp is the Sweet-Parker reconnection rate (Eq. ([!])); 
and Vcsp is the corrected Sweet-Parker reconnection speed given by Eq. fl5"5|). V rec for H 2 is 
placed in brackets as it allows for heating of the reconnection layer. 
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Fig. 1. — A schematic of a reconnection region. Magnetic reconnection is occuring in a layer 
of width A. Gas velocities are marked with thick lines for ions and thin lines for neutrals. 
The gas outflow is confined to a region of width L n . The x axis is horizontal and the y axis 
is vertical. 



